Heterogeneous associations between interleukin-6 receptor variants and phenotypes across ancestries and implications for therapy.

The Phenome-Wide Association Study (PheWAS) is increasingly used to broadly screen for potential treatment effects, e.g., IL6R variant as a proxy for IL6R antagonists. This approach offers an opportunity to address the limited power in clinical trials to study differential treatment effects across patient subgroups. However, limited methods exist to efficiently test for differences across subgroups in the thousands of multiple comparisons generated as part of a PheWAS. In this study, we developed an approach that maximizes the power to test for heterogeneous genotype-phenotype associations and applied this approach to an IL6R PheWAS among individuals of African (AFR) and European (EUR) ancestries. We identified 29 traits with differences in IL6R variant-phenotype associations, including a lower risk of type 2 diabetes in AFR (OR 0.96) vs EUR (OR 1.0, p-value for heterogeneity = 8.5 × 10-3), and higher white blood cell count (p-value for heterogeneity = 8.5 × 10-131). These data suggest a more salutary effect of IL6R blockade for T2D among individuals of AFR vs EUR ancestry and provide data to inform ongoing clinical trials targeting IL6 for an expanding number of conditions. Moreover, the method to test for heterogeneity of associations can be applied broadly to other large-scale genotype-phenotype screens in diverse populations.

Neurological diagnoses in hospitalized COVID-19 patients associated with adverse outcomes: A multinational cohort study.

Few studies examining the patient outcomes of concurrent neurological manifestations during acute COVID-19 leveraged multinational cohorts of adults and children or distinguished between central and peripheral nervous system (CNS vs. PNS) involvement. Using a federated multinational network in which local clinicians and informatics experts curated the electronic health records data, we evaluated the risk of prolonged hospitalization and mortality in hospitalized COVID-19 patients from 21 healthcare systems across 7 countries. For adults, we used a federated learning approach whereby we ran Cox proportional hazard models locally at each healthcare system and performed a meta-analysis on the aggregated results to estimate the overall risk of adverse outcomes across our geographically diverse populations. For children, we reported descriptive statistics separately due to their low frequency of neurological involvement and poor outcomes. Among the 106,229 hospitalized COVID-19 patients (104,031 patients ≥18 years; 2,198 patients <18 years, January 2020-October 2021), 15,101 (14%) had at least one CNS diagnosis, while 2,788 (3%) had at least one PNS diagnosis. After controlling for demographics and pre-existing conditions, adults with CNS involvement had longer hospital stay (11 versus 6 days) and greater risk of (Hazard Ratio = 1.78) and faster time to death (12 versus 24 days) than patients with no neurological condition (NNC) during acute COVID-19 hospitalization. Adults with PNS involvement also had longer hospital stay but lower risk of mortality than the NNC group. Although children had a low frequency of neurological involvement during COVID-19 hospitalization, a substantially higher proportion of children with CNS involvement died compared to those with NNC (6% vs 1%). Overall, patients with concurrent CNS manifestation during acute COVID-19 hospitalization faced greater risks for adverse clinical outcomes than patients without any neurological diagnosis. Our global informatics framework using a federated approach (versus a centralized data collection approach) has utility for clinical discovery beyond COVID-19.

Type 1 Diabetes Genetic Risk in 109,954 Veterans With Adult-Onset Diabetes: The Million Veteran Program (MVP).

OBJECTIVE: To characterize high type 1 diabetes (T1D) genetic risk in a population where type 2 diabetes (T2D) predominates. RESEARCH DESIGN AND METHODS: Characteristics typically associated with T1D were assessed in 109,594 Million Veteran Program participants with adult-onset diabetes, 2011-2021, who had T1D genetic risk scores (GRS) defined as low (0 to <45%), medium (45 to <90%), high (90 to <95%), or highest (≥95%). RESULTS: T1D characteristics increased progressively with higher genetic risk (P < 0.001 for trend). A GRS ≥ 90% was more common with diabetes diagnoses before age 40 years, but 95% of those participants were diagnosed at age ≥40 years, and they resembled T2D in mean age (64.3 years) and BMI (32.3 kg/m2). Compared with the low risk group, the highest-risk group was more likely to have diabetic ketoacidosis (low 0.9% vs. highest GRS 3.7%), hypoglycemia prompting emergency visits (3.7% vs. 5.8%), outpatient plasma glucose <50 mg/dL (7.5% vs. 13.4%), a shorter median time to start insulin (3.5 vs. 1.4 years), use of a T1D diagnostic code (16.3% vs. 28.1%), low C-peptide levels if tested (1.8% vs. 32.4%), and glutamic acid decarboxylase antibodies (6.9% vs. 45.2%), all P < 0.001. CONCLUSIONS: Characteristics associated with T1D were increased with higher genetic risk, and especially with the top 10% of risk. However, the age and BMI of those participants resemble people with T2D, and a substantial proportion did not have diagnostic testing or use of T1D diagnostic codes. T1D genetic screening could be used to aid identification of adult-onset T1D in settings in which T2D predominates.

Long-Term Body Mass Index Variability and Adverse Cardiovascular Outcomes.

Importance: Body mass index (BMI; calculated as weight in kilograms divided by height in meters squared) is a commonly used estimate of obesity, which is a complex trait affected by genetic and lifestyle factors. Marked weight gain and loss could be associated with adverse biological processes. Objective: To evaluate the association between BMI variability and incident cardiovascular disease (CVD) events in 2 distinct cohorts. Design, Setting, and Participants: This cohort study used data from the Million Veteran Program (MVP) between 2011 and 2018 and participants in the UK Biobank (UKB) enrolled between 2006 and 2010. Participants were followed up for a median of 3.8 (5th-95th percentile, 3.5) years. Participants with baseline CVD or cancer were excluded. Data were analyzed from September 2022 and September 2023. Exposure: BMI variability was calculated by the retrospective SD and coefficient of variation (CV) using multiple clinical BMI measurements up to the baseline. Main Outcomes and Measures: The main outcome was incident composite CVD events (incident nonfatal myocardial infarction, acute ischemic stroke, and cardiovascular death), assessed using Cox proportional hazards modeling after adjustment for CVD risk factors, including age, sex, mean BMI, systolic blood pressure, total cholesterol, high-density lipoprotein cholesterol, smoking status, diabetes status, and statin use. Secondary analysis assessed whether associations were dependent on the polygenic score of BMI. Results: Among 92 363 US veterans in the MVP cohort (81 675 [88%] male; mean [SD] age, 56.7 [14.1] years), there were 9695 Hispanic participants, 22 488 non-Hispanic Black participants, and 60 180 non-Hispanic White participants. A total of 4811 composite CVD events were observed from 2011 to 2018. The CV of BMI was associated with 16% higher risk for composite CVD across all groups (hazard ratio [HR], 1.16; 95% CI, 1.13-1.19). These associations were unchanged among subgroups and after adjustment for the polygenic score of BMI. The UKB cohort included 65 047 individuals (mean [SD] age, 57.30 (7.77) years; 38 065 [59%] female) and had 6934 composite CVD events. Each 1-SD increase in BMI variability in the UKB cohort was associated with 8% increased risk of cardiovascular death (HR, 1.08; 95% CI, 1.04-1.11). Conclusions and Relevance: This cohort study found that among US veterans, higher BMI variability was a significant risk marker associated with adverse cardiovascular events independent of mean BMI across major racial and ethnic groups. Results were consistent in the UKB for the cardiovascular death end point. Further studies should investigate the phenotype of high BMI variability.

Centralized Interactive Phenomics Resource: an integrated online phenomics knowledgebase for health data users.

OBJECTIVE: Development of clinical phenotypes from electronic health records (EHRs) can be resource intensive. Several phenotype libraries have been created to facilitate reuse of definitions. However, these platforms vary in target audience and utility. We describe the development of the Centralized Interactive Phenomics Resource (CIPHER) knowledgebase, a comprehensive public-facing phenotype library, which aims to facilitate clinical and health services research. MATERIALS AND METHODS: The platform was designed to collect and catalog EHR-based computable phenotype algorithms from any healthcare system, scale metadata management, facilitate phenotype discovery, and allow for integration of tools and user workflows. Phenomics experts were engaged in the development and testing of the site. RESULTS: The knowledgebase stores phenotype metadata using the CIPHER standard, and definitions are accessible through complex searching. Phenotypes are contributed to the knowledgebase via webform, allowing metadata validation. Data visualization tools linking to the knowledgebase enhance user interaction with content and accelerate phenotype development. DISCUSSION: The CIPHER knowledgebase was developed in the largest healthcare system in the United States and piloted with external partners. The design of the CIPHER website supports a variety of front-end tools and features to facilitate phenotype development and reuse. Health data users are encouraged to contribute their algorithms to the knowledgebase for wider dissemination to the research community, and to use the platform as a springboard for phenotyping. CONCLUSION: CIPHER is a public resource for all health data users available at https://phenomics.va.ornl.gov/ which facilitates phenotype reuse, development, and dissemination of phenotyping knowledge.

Lipid levels and risk of acute pancreatitis using bidirectional Mendelian randomization.

Previous studies found lipid levels, especially triglycerides (TG), are associated with acute pancreatitis, but their causalities and bi-directions were not fully examined. We determined whether abnormal levels of TG, high-density lipoprotein cholesterol (HDL-C), and low-density lipoprotein cholesterol (LDL-C) are precursors and/or consequences of acute pancreatitis using bidirectional two-sample Mendelian randomization (MR) with two non-overlapping genome-wide association study (GWAS) summary statistics for lipid levels and acute pancreatitis. We found phenotypic associations that both higher TG levels and lower HDL-C levels contributed to increased risk of acute pancreatitis. Our GWAS meta-analysis of acute pancreatitis identified seven independent signals. Genetically predicted TG was positively associated with acute pancreatitis when using the variants specifically associated with TG using univariable MR [Odds ratio (OR), 95% CI 2.02, 1.22-3.31], but the reversed direction from acute pancreatitis to TG was not observed (mean difference = 0.003, SE = 0.002, P-value = 0.138). However, a bidirectional relationship of HDL-C and acute pancreatitis was observed: A 1-SD increment of genetically predicted HDL-C was associated with lower risk of acute pancreatitis (OR, 95% CI 0.84, 0.76-0.92) and genetically predisposed individuals with acute pancreatitis have, on average, 0.005 SD lower HDL-C (mean difference = - 0.005, SE = 0.002, P-value = 0.004). Our MR analysis confirms the evidence of TG as a risk factor of acute pancreatitis but not a consequence. A potential bidirectional relationship of HDL-C and acute pancreatitis occurs and raises the prospect of HDL-C modulation in the acute pancreatitis prevention and treatment.

Novel Polygenic Risk Score and Established Clinical Risk Factors for Risk Estimation of Aortic Stenosis.

Importance: Polygenic risk scores (PRSs) have proven to be as strong as or stronger than established clinical risk factors for many cardiovascular phenotypes. Whether this is true for aortic stenosis remains unknown. Objective: To develop a novel aortic stenosis PRS and compare its aortic stenosis risk estimation to established clinical risk factors. Design, Setting, and Participants: This was a longitudinal cohort study using data from the Million Veteran Program (MVP; 2011-2020), UK Biobank (2006-2010), and 6 Thrombolysis in Myocardial Infarction (TIMI) trials, including DECLARE-TIMI 58 (2013-2018), FOURIER (TIMI 59; 2013-2017), PEGASUS-TIMI 54 (2010-2014), SAVOR-TIMI 53 (2010-2013), SOLID-TIMI 52 (2009-2014), and ENGAGE AF-TIMI 48 (2008-2013), which were a mix of population-based and randomized clinical trials. Individuals from UK Biobank and the MVP meeting a previously validated case/control definition for aortic stenosis were included. All individuals from TIMI trials were included unless they had a documented preexisting aortic valve replacement. Analysis took place from January 2022 to December 2023. Exposures: PRS for aortic stenosis (developed using data from MVP and validated in UK Biobank) and other previously validated cardiovascular PRSs, defined either as a continuous variable or as low (bottom 20%), intermediate, and high (top 20%), and clinical risk factors. Main Outcomes: Aortic stenosis (defined using International Classification of Diseases or Current Procedural Terminology codes in UK Biobank and MVP or safety event data in the TIMI trials). Results: The median (IQR) age in MVP was 67 (57-73) years, and 135 140 of 147 104 participants (92%) were male. The median (IQR) age in the TIMI trials was 66 (54-78) years, and 45 524 of 59 866 participants (71%) were male. The best aortic stenosis PRS incorporated 5 170 041 single-nucleotide variants and was associated with aortic stenosis in both the MVP testing sample (odds ratio, 1.41; 95% CI, 1.37-1.45 per 1 SD PRS; P = 4.6 × 10-116) and TIMI trials (hazard ratio, 1.44; 95% CI, 1.27-1.62 per 1 SD PRS; P = 3.2 × 10-9). Among genetic and clinical risk factors, the aortic stenosis PRS performed comparably to most risk factors besides age, and within a given age range, the combination of clinical and genetic risk factors was additive, providing a 3- to 4-fold increased gradient of risk of aortic stenosis. However, the addition of the aortic stenosis PRS to a model including clinical risk factors only improved risk discrimination of aortic stenosis by 0.01 to 0.02 (C index in MVP: 0.78 with clinical risk factors, 0.79 with risk factors and aortic stenosis PRS; C index in TIMI: 0.71 with clinical risk factors, 0.73 with risk factors and aortic stenosis PRS). Conclusions: This study developed and validated 1 of the first aortic stenosis PRSs. While aortic stenosis genetic risk was independent from clinical risk factors and performed comparably to all other risk factors besides age, genetic risk resulted in only a small improvement in overall aortic stenosis risk discrimination beyond age and clinical risk factors. This work sets the stage for further development of an aortic stenosis PRS.

Protein identification for stroke progression via Mendelian Randomization in Million Veteran Program and UK Biobank.

Background: Individuals who have experienced a stroke, or transient ischemic attack, face a heightened risk of future cardiovascular events. Identification of genetic and molecular risk factors for subsequent cardiovascular outcomes may identify effective therapeutic targets to improve prognosis after an incident stroke. Methods: We performed genome-wide association studies (GWAS) for subsequent major adverse cardiovascular events (MACE) (Ncases=51,929, Ncntrl=39,980) and subsequent arterial ischemic stroke (AIS) Ncases=45,120, Ncntrl=46,789) after first incident stroke within the Million Veteran Program and UK Biobank. We then used genetic variants associated with proteins (pQTLs) to determine the effect of 1,463 plasma protein abundances on subsequent MACE using Mendelian randomization (MR). Results: Two variants were significantly associated with subsequent cardiovascular events: rs76472767 (OR=0.75, 95% CI = 0.64-0.85, p= 3.69×10-08) with subsequent AIS and rs13294166 (OR=1.52, 95% CI = 1.37-1.67, p=3.77×10-08) with subsequent MACE. Using MR, we identified 2 proteins with an effect on subsequent MACE after a stroke: CCL27 (effect OR= 0.77, 95% CI = 0.66-0.88, adj. p=0.05), and TNFRSF14 (effect OR=1.42, 95% CI = 1.24-1.60, adj. p=0.006). These proteins are not associated with incident AIS and are implicated to have a role in inflammation. Conclusions: We found evidence that two proteins with little effect on incident stroke appear to influence subsequent MACE after incident AIS. These associations suggest that inflammation is a contributing factor to subsequent MACE outcomes after incident AIS and highlights potential novel targets.

Incidence and prognostic significance of newly-diagnosed atrial fibrillation among older U.S. veterans hospitalized with COVID-19.

Most prior studies on the prognostic significance of newly-diagnosed atrial fibrillation (AF) in COVID-19 did not differentiate newly-diagnosed AF from pre-existing AF. To determine the association between newly-diagnosed AF and in-hospital and 30-day mortality among regular users of Veterans Health Administration using data linked to Medicare. We identified Veterans aged ≥ 65 years who were hospitalized for ≥ 24 h with COVID-19 from 06/01/2020 to 1/31/2022 and had ≥ 2 primary care visits within 24 months prior to the index hospitalization. We performed multivariable logistic regression analyses to estimate adjusted risks, risk differences (RD), and odds ratios (OR) for the association between newly-diagnosed AF and the mortality outcomes adjusting for patient demographics, baseline comorbidities, and presence of acute organ dysfunction on admission. Of 23,299 patients in the study cohort, 5.3% had newly-diagnosed AF, and 29.2% had pre-existing AF. In newly-diagnosed AF adjusted in-hospital and 30-day mortality were 16.5% and 22.7%, respectively. Newly-diagnosed AF was associated with increased mortality compared to pre-existing AF (in-hospital: OR 2.02, 95% confidence interval [CI] 1.72-2.37; RD 7.58%, 95% CI 5.54-9.62) (30-day: OR 1.86; 95% CI 1.60-2.16; RD 9.04%, 95% CI 6.61-11.5) or no AF (in-hospital: OR 2.24, 95% CI 1.93-2.60; RD 8.40%, 95% CI 6.44-10.4) (30-day: 2.07, 95% CI 1.80-2.37; RD 10.2%, 95% CI 7.89-12.6). There was a smaller association between pre-existing AF and the mortality outcomes. Newly-diagnosed AF is an important prognostic marker for patients hospitalized with COVID-19. Whether prevention or treatment of AF improves clinical outcomes in these patients remains unknown.

Red Meat Intake and the Risk of Cardiovascular Diseases: A Prospective Cohort Study in the Million Veteran Program.

BACKGROUND: Red meat consumption was associated with an increased risk of cardiovascular disease (CVD) in prospective cohort studies and a profile of biomarkers favoring high CVD risk in short-term controlled trials. However, several recent systematic reviews and meta-analyses concluded with no or weak evidence for limiting red meat intake. OBJECTIVES: To prospectively examine the associations between red meat intake and incident CVD in an ongoing cohort study with diverse socioeconomic and racial or ethnic backgrounds. METHODS: Our study included 148,506 participants [17,804 female (12.0%)] who were free of cancer, diabetes, and CVD at baseline from the Million Veteran Program. A food frequency questionnaire measured red meat intakes at baseline. Nonfatal myocardial infarction and acute ischemic stroke were identified through a high-throughput phenotyping algorithm, and fatal CVD events were identified by searching the National Death Index. RESULTS: Comparing the extreme categories of intake, the multivariate-adjusted relative risks of CVD was 1.18 (95% CI: 1.01, 1.38; P-trend < 0.0001) for total red meat, 1.14 (95% CI: 0.96, 1.36; P-trend = 0.01) for unprocessed red meat, and 1.29 (95% CI: 1.04, 1.60; P-trend = 0.003) for processed red meat. We observed a more pronounced positive association between red meat intake and CVD in African American participants than in White participants (P-interaction = 0.01). Replacing 0.5 servings/d of red meat with 0.5 servings/d of nuts, whole grains, and skimmed milk was associated with 14% (RR: 0.86; 95% CI: 0.83, 0.90), 7% (RR: 0.93; 95% CI: 0.89, 0.96), and 4% (RR: 0.96; 95% CI: 0.94, 0.99) lower risks of CVD, respectively. CONCLUSIONS: Red meat consumption is associated with an increased risk of CVD. Our findings support lowering red meat intake and replacing red meat with plant-based protein sources or low-fat dairy foods as a key dietary recommendation for the prevention of CVD.

The Million Veteran Program 1990-1991 Gulf War Era Survey: An Evaluation of Veteran Response, Characteristics, and Representativeness of the Gulf War Era Veteran Population.

To address gaps in understanding the pathophysiology of Gulf War Illness (GWI), the VA Million Veteran Program (MVP) developed and implemented a survey to MVP enrollees who served in the U.S. military during the 1990-1991 Persian Gulf War (GW). Eligible Veterans were invited via mail to complete a survey assessing health conditions as well as GW-specific deployment characteristics and exposures. We evaluated the representativeness of this GW-era cohort relative to the broader population by comparing demographic, military, and health characteristics between respondents and non-respondents, as well as with all GW-era Veterans who have used Veterans Health Administration (VHA) services and the full population of U.S. GW-deployed Veterans. A total of 109,976 MVP GW-era Veterans were invited to participate and 45,270 (41%) returned a completed survey. Respondents were 84% male, 72% White, 8% Hispanic, with a mean age of 61.6 years (SD = 8.5). Respondents were more likely to be older, White, married, better educated, slightly healthier, and have higher socioeconomic status than non-respondents, but reported similar medical conditions and comparable health status. Although generally similar to all GW-era Veterans using VHA services and the full population of U.S. GW Veterans, respondents included higher proportions of women and military officers, and were slightly older. In conclusion, sample characteristics of the MVP GW-era cohort can be considered generally representative of the broader GW-era Veteran population. The sample represents the largest research cohort of GW-era Veterans established to date and provides a uniquely valuable resource for conducting in-depth studies to evaluate health conditions affecting 1990-1991 GW-era Veterans.

LATTE: Label-efficient incident phenotyping from longitudinal electronic health records.

Electronic health record (EHR) data are increasingly used to support real-world evidence studies but are limited by the lack of precise timings of clinical events. Here, we propose a label-efficient incident phenotyping (LATTE) algorithm to accurately annotate the timing of clinical events from longitudinal EHR data. By leveraging the pre-trained semantic embeddings, LATTE selects predictive features and compresses their information into longitudinal visit embeddings through visit attention learning. LATTE models the sequential dependency between the target event and visit embeddings to derive the timings. To improve label efficiency, LATTE constructs longitudinal silver-standard labels from unlabeled patients to perform semi-supervised training. LATTE is evaluated on the onset of type 2 diabetes, heart failure, and relapses of multiple sclerosis. LATTE consistently achieves substantial improvements over benchmark methods while providing high prediction interpretability. The event timings are shown to help discover risk factors of heart failure among patients with rheumatoid arthritis.

High dimensional predictions of suicide risk in 4.2 million US Veterans using ensemble transfer learning.

We present an ensemble transfer learning method to predict suicide from Veterans Affairs (VA) electronic medical records (EMR). A diverse set of base models was trained to predict a binary outcome constructed from reported suicide, suicide attempt, and overdose diagnoses with varying choices of study design and prediction methodology. Each model used twenty cross-sectional and 190 longitudinal variables observed in eight time intervals covering 7.5 years prior to the time of prediction. Ensembles of seven base models were created and fine-tuned with ten variables expected to change with study design and outcome definition in order to predict suicide and combined outcome in a prospective cohort. The ensemble models achieved c-statistics of 0.73 on 2-year suicide risk and 0.83 on the combined outcome when predicting on a prospective cohort of [Formula: see text] 4.2 M veterans. The ensembles rely on nonlinear base models trained using a matched retrospective nested case-control (Rcc) study cohort and show good calibration across a diversity of subgroups, including risk strata, age, sex, race, and level of healthcare utilization. In addition, a linear Rcc base model provided a rich set of biological predictors, including indicators of suicide, substance use disorder, mental health diagnoses and treatments, hypoxia and vascular damage, and demographics.

New statin use, mortality, and first cardiovascular events in older US Veterans by frailty status.

BACKGROUND: Statins are part of long-term medical regimens for many older adults. Whether frailty modifies the protective relationship between statins, mortality, and major adverse cardiovascular events (MACE) is unknown. METHODS: This was a retrospective study of US Veterans ≥65, without CVD or prior statin use seen in 2002-2012, followed through 2017. A 31-item frailty index was used. The co-primary endpoint was all-cause mortality or MACE (MI, stroke/TIA, revascularization, or cardiovascular death). Cox proportional hazards models were developed to evaluate the association of statin use with outcomes; propensity score overlap weighting accounted for confounding by indication. RESULTS: We identified 710,313 Veterans (mean age (SD) 75.3(6.5), 98% male, 89% white); 86,327 (12.1%) were frail. Over mean follow-up of 8 (5) years, there were 48.6 and 72.6 deaths per 1000 person-years (PY) among non-frail statin-users vs nonusers (weighted Incidence Rate Difference (wIRD)/1000 person years (PY), -24.0[95% CI, -24.5 to -23.6]), and 90.4 and 130.4 deaths per 1000PY among frail statin-users vs nonusers (wIRD/1000PY, -40.0[95% CI, -41.8 to -38.2]). There were 51.7 and 60.8 MACE per 1000PY among non-frail statin-users vs nonusers (wIRD/1000PY, -9.1[95% CI, -9.7 to -8.5]), and 88.2 and 102.0 MACE per 1000PY among frail statin-users vs nonusers (wIRD/1000PY, -13.8[95% CI, -16.2 to -11.4]). There were no significant interactions by frailty for statin users vs non-users by either mortality or MACE outcomes, p-interaction 0.770 and 0.319, respectively. Statin use was associated with lower risk of all-cause mortality (HR, 0.61 (0.60-0.61)) and MACE (HR 0.86 (0.85-0.87)). CONCLUSIONS: New statin use is associated with a lower risk of mortality and MACE, independent of frailty. These findings should be confirmed in a randomized clinical trial.

Impact of 8 lifestyle factors on mortality and life expectancy among United States veterans: The Million Veteran Program.

BACKGROUND: Lifestyle medicine has been proposed as a way to address the root causes of chronic disease and their associated health care costs. OBJECTIVE: This study aimed to estimate mortality risk and longevity associated with individual lifestyle factors and comprehensive lifestyle therapy. METHODS: Age- and sex-specific mortality rates were calculated on the basis of 719,147 veterans aged 40-99 y enrolled in the Veteran Affairs Million Veteran Program (2011-2019). Hazard ratios and estimated increase in life expectancy were examined among a subgroup of 276,132 veterans with complete data on 8 lifestyle factors at baseline. The 8 lifestyle factors included never smoking, physical activity, no excessive alcohol consumption, restorative sleep, nutrition, stress management, social connections, and no opioid use disorder. RESULTS: On the basis of 1.12 million person-years of follow-up, 34,247 deaths were recorded. Among veterans who adopted 1, 2, 3, 4, 5, 6, 7, and 8 lifestyle factors, the adjusted hazard ratios for mortality were 0.74 (0.60-0.90), 0.60 (95% CI: 0.49, 0.73), 0.50 (95% CI: 0.41, 0.61), 0.43 (95% CI: 0.35, 0.52), 0.35 (95% CI: 0.29, 0.43), 0.27 (95% CI: 0.22, 0.33), 0.21 (95% CI: 0.17, 0.26), and 0.13 (95% CI: 0.10, 0.16), respectively, as compared with veterans with no adopted lifestyle factors. The estimated life expectancy at age 40 y was 23.0, 26.5, 28.8, 30.8, 32.7, 35.1, 38.3, 41.3, and 47.0 y among males and 27.0, 28.8, 33.1, 38.0, 39.2, 41.4, 43.8, 46.3, and 47.5 y for females who adopted 0, 1, 2, 3, 4, 5, 6, 7, and 8 lifestyle factors, respectively. The difference in life expectancy at age 40 y was 24.0 y for male veterans and 20.5 y for female veterans when comparing adoption of 8-9 lifestyle factors. CONCLUSIONS: A combination of 8 lifestyle factors is associated with a significantly lower risk of premature mortality and an estimated prolonged life expectancy.

Comparison of the Test-negative Design and Cohort Design With Explicit Target Trial Emulation for Evaluating COVID-19 Vaccine Effectiveness.

BACKGROUND: Observational studies are used for estimating vaccine effectiveness under real-world conditions. The practical performance of two common approaches-cohort and test-negative designs-need to be compared for COVID-19 vaccines. METHODS: We compared the cohort and test-negative designs to estimate the effectiveness of the BNT162b2 vaccine against COVID-19 outcomes using nationwide data from the United States Department of Veterans Affairs. Specifically, we (1) explicitly emulated a target trial using follow-up data and evaluated the potential for confounding using negative controls and benchmarking to a randomized trial, (2) performed case-control sampling of the cohort to confirm empirically that the same estimate is obtained, (3) further restricted the sampling to person-days with a test, and (4) implemented additional features of a test-negative design. We also compared their performance in limited datasets. RESULTS: Estimated BNT162b2 vaccine effectiveness was similar under all four designs. Empirical results suggested limited residual confounding by healthcare-seeking behavior. Analyses in limited datasets showed evidence of residual confounding, with estimates biased downward in the cohort design and upward in the test-negative design. CONCLUSION: Vaccine effectiveness estimates under a cohort design with explicit target trial emulation and a test-negative design were similar when using rich information from the VA healthcare system, but diverged in opposite directions when using a limited dataset. In settings like ours with sufficient information on confounders and other key variables, the cohort design with explicit target trial emulation may be preferable as a principled approach that allows estimation of absolute risks and facilitates interpretation of effect estimates.

Rates of ICD-10 Code U09.9 Documentation and Clinical Characteristics of VA Patients With Post-COVID-19 Condition.

Importance: A significant proportion of SARS-CoV-2 infected individuals experience post-COVID-19 condition months after initial infection. Objective: To determine the rates, clinical setting, risk factors, and symptoms associated with the documentation of International Statistical Classification of Diseases Tenth Revision (ICD-10), code U09.9 for post-COVID-19 condition after acute infection. Design, Setting, and Participants: This retrospective cohort study was performed within the US Department of Veterans Affairs (VA) health care system. Veterans with a positive SARS-CoV-2 test result between October 1, 2021, the date ICD-10 code U09.9 was introduced, and January 31, 2023 (n = 388 980), and a randomly selected subsample of patients with the U09.9 code (n = 350) whose symptom prevalence was assessed by systematic medical record review, were included in the analysis. Exposure: Positive SARS-CoV-2 test result. Main Outcomes and Measures: Rates, clinical setting, risk factors, and symptoms associated with ICD-10 code U09.9 in the medical record. Results: Among the 388 980 persons with a positive SARS-CoV-2 test, the mean (SD) age was 61.4 (16.1) years; 87.3% were men. In terms of race and ethnicity, 0.8% were American Indian or Alaska Native, 1.4% were Asian, 20.7% were Black, 9.3% were Hispanic or Latino, 1.0% were Native Hawaiian or Other Pacific Islander; and 67.8% were White. Cumulative incidence of U09.9 documentation was 4.79% (95% CI, 4.73%-4.87%) at 6 months and 5.28% (95% CI, 5.21%-5.36%) at 12 months after infection. Factors independently associated with U09.9 documentation included older age, female sex, Hispanic or Latino ethnicity, comorbidity burden, and severe acute infection manifesting by symptoms, hospitalization, or ventilation. Primary vaccination (adjusted hazard ratio [AHR], 0.80 [95% CI, 0.78-0.83]) and booster vaccination (AHR, 0.66 [95% CI, 0.64-0.69]) were associated with a lower likelihood of U09.9 documentation. Marked differences by geographic region and facility in U09.9 code documentation may reflect local screening and care practices. Among the 350 patients undergoing systematic medical record review, the most common symptoms documented in the medical records among patients with the U09.9 code were shortness of breath (130 [37.1%]), fatigue or exhaustion (78 [22.3%]), cough (63 [18.0%]), reduced cognitive function or brain fog (22 [6.3%]), and change in smell and/or taste (20 [5.7%]). Conclusions and Relevance: In this cohort study of 388 980 veterans, documentation of ICD-10 code U09.9 had marked regional and facility-level variability. Strong risk factors for U09.9 documentation were identified, while vaccination appeared to be protective. Accurate and consistent documentation of U09.9 is needed to maximize its utility in tracking patients for clinical care and research. Future studies should examine the long-term trajectory of individuals with U09.9 documentation.

Moderate alcohol consumption on the risk of stroke in the Million Veteran Program.

BACKGROUND: There is inconsistent evidence on the association of moderate alcohol consumption and stroke risk in the general population and is not well studied among U.S. Veterans. Furthermore, it is unclear whether primarily drinking beer, wine, or liquor is associated with a difference in stroke risk. METHODS: The study included 185,323 Million Veteran Program participants who self-reported alcohol consumption on the Lifestyle Survey. Moderate consumption was defined as 1-2 drinks/day and beverage preference of beer, wine or liquor was defined if ≥ 50% of total drinks consumed were from a single type of beverage. Strokes were defined using ICD-9 and ICD-10 codes from the participants' electronic health record. RESULTS: The mean (sd) age of the sample was 64 (13) years and 11% were women. We observed 4,339 (94% ischemic; 6% hemorrhagic) strokes over a median follow-up of 5.2 years. In Cox models adjusted for age, sex, race, education, income, body mass index, smoking, exercise, diet, cholesterol, prevalent diabetes, prevalent hypertension, lipid-lowering medication, antihypertensive medication, and diabetes medication, moderate alcohol consumption (1-2 drinks/day) was associated with a 22% lower risk of total stroke compared with never drinking [Hazards ratio (HR) 95% confidence interval (CI): 0.78 (0.67, 0.92)]. When stratifying by stroke type, we observed a similar protective association with moderate consumption and ischemic stroke [HR (95% CI): 0.76 (0.65, 0.90)], but a non-statistically significant higher risk of hemorrhagic stroke [HR (95% CI): 1.29 (0.64, 2.61)]. We did not observe a difference in ischemic or hemorrhagic stroke risk among those who preferred beer, liquor or wine vs. no beverage preference. When stratifying by prior number of hospital visits (≤ 15, 16-33, 34-64, ≥ 65) as a proxy for health status, we observed attenuation of the protective association with greater number of visits [HR (95% CI): 0.87 (0.63, 1.19) for ≥ 65 visits vs. 0.80 (0.59, 1.08) for ≤ 15 visits]. CONCLUSIONS: We observed a lower risk of ischemic stroke, but not hemorrhagic stroke with moderate alcohol consumption and did not observe substantial differences in risk by beverage preference among a sample of U.S. Veterans. Healthy user bias of moderate alcohol consumption may be driving some of the observed protective association.

Characterizing prostate cancer risk through multi-ancestry genome-wide discovery of 187 novel risk variants.

The transferability and clinical value of genetic risk scores (GRSs) across populations remain limited due to an imbalance in genetic studies across ancestrally diverse populations. Here we conducted a multi-ancestry genome-wide association study of 156,319 prostate cancer cases and 788,443 controls of European, African, Asian and Hispanic men, reflecting a 57% increase in the number of non-European cases over previous prostate cancer genome-wide association studies. We identified 187 novel risk variants for prostate cancer, increasing the total number of risk variants to 451. An externally replicated multi-ancestry GRS was associated with risk that ranged from 1.8 (per standard deviation) in African ancestry men to 2.2 in European ancestry men. The GRS was associated with a greater risk of aggressive versus non-aggressive disease in men of African ancestry (P = 0.03). Our study presents novel prostate cancer susceptibility loci and a GRS with effective risk stratification across ancestry groups.

Knowledge-Driven Online Multimodal Automated Phenotyping System.

Though electronic health record (EHR) systems are a rich repository of clinical information with large potential, the use of EHR-based phenotyping algorithms is often hindered by inaccurate diagnostic records, the presence of many irrelevant features, and the requirement for a human-labeled training set. In this paper, we describe a knowledge-driven online multimodal automated phenotyping (KOMAP) system that i) generates a list of informative features by an online narrative and codified feature search engine (ONCE) and ii) enables the training of a multimodal phenotyping algorithm based on summary data. Powered by composite knowledge from multiple EHR sources, online article corpora, and a large language model, features selected by ONCE show high concordance with the state-of-the-art AI models (GPT4 and ChatGPT) and encourage large-scale phenotyping by providing a smaller but highly relevant feature set. Validation of the KOMAP system across four healthcare centers suggests that it can generate efficient phenotyping algorithms with robust performance. Compared to other methods requiring patient-level inputs and gold-standard labels, the fully online KOMAP provides a significant opportunity to enable multi-center collaboration.